Study of the pp — > ppn + n Reaction in the Low-Energy Tail of the Roper Resonance 
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Exclusive measurements of the pp — > pp-K + -R~ reaction have been carried out at T p — 775 MeV 
at CELSIUS using the PROMICE/WASA setup. Together with data obtained at lower energy they 
point to a dominance of the Roper excitation in this process. From the observed interference of its 
decay routes N* — > Na and N* — > An — > Na their energy-dependent relative branching ratio is 
determined. 

PACS numbers: 13.75.-n, 14.20.Gk, 25.40.Ve 



The Roper resonance A* (1440) with I{J P ) = 
1/2(1/2+) is presently known as the second excited state 
of the nucleon pj. In contrast to the first excited 
state, the A(1232), and also higher-lying resonances the 
A* (1440) is still poorly understood both theoretically 
and experimentally. Since it is hardly excited by elec- 
tromagnetic probes and has the same quantum numbers 
as the nucleon, it has been interpreted as the breath- 
ing mode monopole excitation of the nucleon. A recent 
theoretical work 2| Q finds the Roper excitation to rest 
solely on meson- nucleon dynamics, whereas another re- 
cent investigation £| proposes it to be actually two reso- 
nances with one being the breathing mode and the other 
one a A excitation built on top of the A(1232). In all 
these aspects the decay modes of the Roper resonance 
into the Nirir channel play a crucial role. The simplest 
decay is A* — > N(irn)i = i = o := Na, i.e., the decay into 
the a channel. A competitive and according to present 
knowledge Q actually much stronger decay channel is 
the sequential Roper decay via the A(1232) resonance, 
A* — > A-7T. However, this decay channel is not very well 
defined, in particular not orthogonal to the Act channel, 
since the A is also unstable and decays nearly as fast as 
the Roper does. In fact, most of this decay will end up 
again in the Na channel and thus will interfere with the 
direct N* — > Na decay. 



Valencia group exhibits this reaction to be unique in 
the sense that it selectively provides the excitation mode 
"ct" N — » A* (where "a" stands for the a exchange), 
which is not accessible in any other basic reaction process 
leading to the Roper excitation. 

In this work we present new data from an exclusive 



measurement of the pp 



ppTT + Tr 



reaction at T B = 775 



In a previous work, the first exclusive measurement of 
the pp — * pptt + tt~ reaction at T p — 750 MeV we 
have shown that at energies not far above threshold this 
reaction can be well described by dominant a exchange 
in the initial N N collision with subsequent excitation of 
the Roper resonance in one of the nucleons. This result, 
which is in agreement with theoretical predictions of the 



MeV. Together with the data at T p = 750 MeV they are 
analysed with particular emphasis on the interference of 
the decay routes A* -> Act and A* -» Att -> Na. The 
data have been taken at the CELSIUS storage ring us- 
ing the PROMICE/WASA detector with a cluster jet H 2 
target ,7] . Protons and ir + particles have been registered 
in the forward detector part covering the polar angles 
4° < 9lq6 < 21°. The particles have been identified by 
the AE — E method, the stopped tt + particles in addi- 
tion by their delayed pulse from subsequent muon decay. 
From the measured four-momenta of the two registered 

protons and the identified 7r + the full ppir + Tr~ events 
have been reconstructed by kinematical fits with one 
overconstraint. Detector efficiencies and acceptance have 
been obtained from Monte-Carlo (MC) simulations of the 
detector response Q • The absolute normalization of the 
data has been obtained by monitoring the luminosity of 
the experiment by the simultaneous measurement of the 
elastic scattering and its comparison to data from liter- 
ature 0- The total cross section for the pp — * pp7r + ir~ 
reaction obtained for T p = 775 MeV is a to t — 2.2(5)^ib 
and shown in Fig. 1, together with previous results 
HQS mill El 111 El Our value is an order of magni- 
tude below the bubble chamber results [13J , in agreement 
with our findings at lower energies jm . The estimated 
uncertainty of about 20% is due to y| statistical uncer- 
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FIG. 2: Definition of the different scattering angles in the sub- 
system of particles, here for the case of two particles (pl,p2) 
resulting from the reaction in the overall centre of mass sys- 
tem. For simplicity, the figure shows a non-relativistic con- 
struction. The other particles are not shown. For details see 
text. 



FIG. 1: Energy excitation function of the integral cross sec- 
tion of the pp -» pptt+tt" reaction. The WASA/PROMICE 
data (black points) from this work and Ref. [5] are compared 
with previous data (open symbols) [10-15] and predictions [6] 
with two different parameter sets (solid and dashed line), with 
and without pp final state interaction (upper and lower lines 
respectively) and a phase space distribution adjusted to our 
data (dotted line). 



taintics in the collected pp elastic (1%) and ppn + ir~ (5%) 
events as well as systematic uncertainties in the selection 
of elastic (8%) andpp7r + 7r~ (12%) events, uncertainty 
in the lifetime of the data acquisition system (8%) and 
uncertainties in the extrapolation to full solid angle (6%). 
For T p = 750 MeV we show two values for av ot . The up- 
per one is the previously published value |5j, the lower 
one has been derived from a subsample of those data us- 
ing the same event selection criteria as applied now for 
the 775 MeV data. In order to test the robustness of the 
analysis the event selection criteria have been slightly 
modified compared to the ones used previously Q . How- 
ever, within uncertainties both values agree with each 
other. 

Differential cross sections are shown in Figs. 3-5, which 
will be discussed in the following. As pointed out in 
Ref. 5], the proton angular distribution in the over- 
all center of mass system (cms) is governed by the me- 
son exchange between the colliding protons. In partic- 
ular, it should be described by a-exchange leading to 
cr(Q p ) ~ 1 — acos 2 (8 p ) with a > given by the am- 
plitude for a exchange while a < would be typical for 
7r-exchange. The data at T p — 750 MeV were well de- 
scribed by this ansatz. This holds also for the new data 
at T p = 775 MeV @ . Instead of looking at the angular 
distribution in the overall cms it appears more instructive 
to look at the angular distributions in the pp subsystem. 
The definition of angles is illustrated in Fig. 2. Let us 
denote the scattering angles of particles 1 and 2 in the 
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FIG. 3: Proton angular distributions in the pp subsystem at 
T p = 775 MeV. In the left diagram the proton emission angle 
<dp P is taken relative to the protons' summary momentum in 
the overall cms. Shaded areas give the phase space distri- 
bution, whereas dashed and solid lines show MC simulations 
according to ansatz (1) and (2), respectively. 



overall cms by Pl and Q P2 , respectively, and the scat- 
tering angle of the center of mass motion of both protons 
(summed momenta) in the overall cms by Q PlP2 ■ Within 
the rest frame of the two particles (p\P2 subsystem) then 
two angles can be defined, the scattering angle of P2 ei- 
ther with respect to the beam axis, Op^ 2 , or with respect 
to the summed momenta of p\ and pi in the overall cms. 
The latter angle is denoted by 9p^ P2 . Note that due to 
the indistinguishability of the two protons all these angu- 
lar distributions have to be symmetric about 90°. Since 
in this reaction the two protons are emitted dominantly 
back-to-back in the overall cms (see distributions of the 
opening angle S pp between the two protons in Refs. 
the <d pp distribution (Fig. 3) is very close to the P dis- 
tribution and 

hence also exhibits a (1 — a cos 2 O pp ) dependence. Shown 

in Fig. 3 is in addition the Q pp distribution. This dis- 
tribution reflects the scattering situation of the two pro- 
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tons within their subsystem. The observed distribution 
is isotropic, i.e. we indeed see the outgoing protons to be 
in relative s-wave, as anticipated in Ref . |5j ■ 

We note in passing that also the pion angular distri- 
bution in the overall cms |8j is flat as well — as was the 
case at T p — 750 MeV, too. Since here only the pro- 
cess NN — > AA — > NNtttt leads to nonisotropic angu- 
lar distributions [EIEQj we conclude that this process, 
which is expected to contribute substantially at energies 
T p > 1000 MeV 6] is not yet of relevance at T p = 775 
MeV. 

To see whether the reaction proceeds via N* excita- 
tion, we inspect the measured distribution of the pn + n~ 
invariant mass M pn + n - (Fig. 4). Compared to phase 
space the data are substantially enhanced near the high- 
energy end, compatible with the low-energy tail of the 
N* excitation and reproduced by the appropriate cal- 
culations for N* excitation, as will be discussed in de- 
tail below. As in Ref. 5] we conclude that the process 
u cr"N — > N* is indeed the one which drives the reaction 
pp — > ppir + ir~ at the energies considered here. 

We now examine the N* decay process as exhibited 
by the data. In the analysis of the 750 MeV data two 
versions for the N* decay amplitude have been proposed 



T p = 750 MeV 



T p = 775 MeV 



5] 



A ~ 1 + c ki • k 2 (3D A ++ + D A o) 



and 



A~ (1 + c' kx-k 2 )D 



A + + 



(1) 



(2) 



In the full reaction amplitude this factor A complements 
the propagators for a exchange and N* excitation as 
well as the expression describing the final state inter- 
action between the two outgoing protons in relative s- 
wave. Here D A ++ = l/(M p „+ - M A ++ + ^r A++ ) and 
_Da° defined analogously, are the A propagators. The 
constant 1 stands for the process N* — > Na and the 
scalar product ki ■ k 2 of the pion momenta ki and k 2 
for the double p-w&ve decay of the route N* — > An — » 
Na. For simplicity we have neglected the spinflip term 
is ■ (ki x k 2 )(3DA++ — D A o) in this decay channel with 
s being the nucleon spin. This term describes the transi- 
tion N* — > Air — ^ N(tvtt)i—i—i. It does not interfere with 
the other terms, is smaller than the ki • k 2 term by a fac- 
tor of 16 in the cross section and hence does not influence 
significantly the conclusions in this paper. Whereas in 
Refs. Hi' this scalar product has been calculated in the 
overall cms, we here use the more appropriate N* sys- 
tem. We note, however, that the difference is tiny, since 
near threshold the static limit approximation is reason- 
ably valid. 

Ansatz (1) represents the leading term of the two-pion 
decay of the Roper resonance as worked out by the Va- 
lencia group 0. The constant c (and correspondingly 
d in ansatz (2)) gives the relative strength between the 
two decay routes and is treated in the following as the 
parameter to be adjusted to the data, which will enable 
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FIG. 4: Invariant masses M pn + n - and M v + W - as well as the 
opening angle 5^+^- between the two pions for both 750 (left) 
and 775 MeV (right) beam energies. The data are shown in 
comparison to phase space (shaded area) and MC simulations 
for pure decays N* — > Na (dotted), TV* — > An (dashed dot- 
ted) and their interference with c' = —37 (dashed), and -61 
(solid) (GeV/c)" 2 using eq. (2). 



us to deduce the relative branching ratio of the two decay 
routes in question. With c being adjusted appropriately 
we get a quantitative description of the data both for 
T p = 750 MeV @ and 775 MeV (see Figs. 3-5) with the 
exception of the distributions for the invariant masses 
A/ p7r +,M p7r - and 0£+' r (dashed lines in Fig. 5). As 
shown in Ref. [j| ansatz (2) is able to heal this deficiency 
(solid lines in Fig. 5) without destroying the good agree- 
ment in the other observables. We admit, however, that 
eq. (2) having the A propagator multiplying also the con- 
stant term, is a purely phcnomcnological ad hoc ansatz. 
Apparently it is successful, but its physics contents is not 
(yet) fully understood. As suggested in Ref. it pos- 
sibly accounts effectively for some final state interaction 
effect. Whereas in ansatz (1) the parameter c is defined 
as energy-independent, since all dynamics is taken into 
account explicitly, the situation is not so clear with our 
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phenomenological ansatz (2); if some final-state interac- 
tion is absorbed here, an energy dependence of c' cannot 
be excluded a priori. 

In Fig. 4 we compare calculations assuming different 
mixing scenarios for the two N* decay routes to the data 
for those observables which are most sensitive to the N* 
decays. Shown are the distributions of invariant masses 
Mp^+K- and M ff + ff - as well as 5^+^- = 

(ki,k 2 ), i.e. the opening angle between the two 
pion momenta in the overall cms. The latter distribu- 
tion directly reflects the squared decay amplitudes (1) 
and (2), respectively, averaged over all possible mod- 
uli of pion momenta at given 5^+ n -, i.e., a(6 7r + 7r -) ~ 



(1 + 6cos<5„. 



with the mixing coefficient b. 



In case of eq. (2) we have b — d < kik% > where the 
brackets denote the average over all possible combina- 
tions. For 1 the distribution cr((5 7r + 7r -) is essentially 
linear in 6, whereas this dependence gets quadratic for 
b 1. Shown in Fig. 4 are calculations for pure phase 
space and for transitions via either the N* — » Na route 
(b = 0) or the N* — > An — > Na route. In order to il- 
lustrate the sensitivity of the data to the mixing of both 
routes, calculations are also shown with c' = —37 and 
—61 (GeV/c) -2 corresponding to b = —0.20 and —0.33, 
respectively, at T p — 750 MeV. The negative sign of the 
coefficients reflects the destructive interference between 
both terms, which is required by the data. If we fit b for 
best reproduction of the data we obtain b = —0.27(2) for 
T p = 750 MeV and b = -0.32(1) for T p = 775 MeV, or 
c' = -50(4) and -53(3) (GeV/c) -2 , respectively. Both 
values agree within their uncertainties, as they should if 
c' is energy independent. I.e. we not only observe the 
proper dependence in the angle <)i(ki, k 2 ), but also in the 
energy implied by fcifc 2 as the beam energy is changed, 
and with it the energy of the N* excitation: for T p = 750 
MeV we have < M N * >= 1264 MeV and for T p = 775 
MeV the average value is < M N , >= 1272 MeV. 

Alternatively, if we use eq. (1) for the description of 
the data we arrive at c = 1.88(8) and 1.96(8) (GeV/c) -1 



for T p = 750 and 775 MeV, respectively. 

Having fitted the parameters c and c', respectively, we 
can determine the ratio of the partial decay widths for the 
routes N* — > Air — » Nirir and N* — > Na in dependence 
of the excited N* mass by 



R(M N ») 
9 , , 



N* 



\M N >) 



i2 J\M A7r \ 2 dM* n+ dM* +w _ 



(3) 



J\M Na \ 2 dM 



,f2 



+ dM 2 + 



with the matrix elements M. Na = 1 and A4a-it — 
ki • k 2 (3£>A++ + -Da ) in case of eq. (1). In case of 
eq. (2) these matrix elements are M n<j — and 
■Mait = ki • k2£>A++- Note that the integral is just 
the integration of the matrix element squared over the 
Dalitz plot in dependence of the invariant masses 
and M 2 + _. The factor 9/8 in eq. (3) is determined by 
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FIG. 5: Invariant masses of pir + and pix~ systems (top and 
middle) as well as scattering angle of the 7r + in the 7T7r-system 
with respect to the summed pion momenta (bottom). The 
left side shows the results for 750 MeV and the right one for 
775 MeV. The shaded areas give phase space distributions, 
whereas dashed and solid lines show MC simulations accord- 
ing to ansatz (1) and (2), respectively. 



isospin coupling coefficients and accounts for the decay 
into channels other than p(7r + 7r~)/ = i = o. If we neglect 
spinflip contributions, then 2/3 of both the p* — > pa de- 
cay and of the p* — > An decay end up in the pn + TT~ 
channel, i.e. the correction factor is unity instead of 9/8 

The results of these calculations are given in Table 1. 
For T p = 750 MeV and T p = 775 MeV both equations 
lead to ratios R(Mn-) which agree within uncertainties. 
In this low-energy tail of the Roper resonance the ra- 
tio turns out to be very small, and is strongly energy- 
dependent as expected from the ki • k 2 dependence of 
Matt- We find #(1272) « 1.5 * i?(1264), i.e. a 50% 
relative increase in the N* — > An route at the higher- 
energy. This increase is essentially due to the increase of 
< fcifo > 2 , which increases by more than 40% by the 25 
MeV increase in the beam energy. The additional energy 
dependence in the N* — > An route due to the A propa- 
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TABLE I: Ratio of the branching ratios for the decays N* — > 
Atv — > Ntttt and N* — * Na in dependence of the excited -/V* 
mass using ansatz (1) and (2), respectively, for the analysis of 
the data at T v = 750 and 775 MeV. The extracted parameters 
c and c' are given in units of (GeV/c)" 1 and (GeV/c) -2 , 
respectively. 

eq- (1) eq. (2) 



(T p = 750 MeV) 



(T p = 775 MeV) 

R(1264) 
(T p = 750 MeV) 

R(1272) 
(T p = 775 MeV) 

R(1371) 
extrapolated 

R(1440) 
extrapolated 

R(1440) 
PDG [1] 



1.96(9) 



-50(4) 
-53(3) 



0.034(4) 0.030(5) 
0.054(6) 0.047(5) 
1.0(1) 0.6(1) 
3.4(3) 1.1(2) 
4(2) 



gator in ansatz (1) is still of minor importance here. In 
case of ansatz (2) such a small energy dependence could 
be compensated easily by a small variation in the pa- 
rameter c' within its statistical uncertainty. Hence with 
regard to this point the two data sets at 750 and 775 MeV 
are not yet able to discriminate between both equations. 



However, the different appearance of the A propaga- 
tor in eqs. (1) and (2) will get discriminative if we go 
to higher energies. To demonstrate this we extrapolate 
R(Mn-) to the nominal resonance pole at 1440 MeV/c 2 
assuming eqs. (1) and (2) to hold also at higher energies 
and taking for c and d the average of the values obtained 
at T p — 750 and 775 MeV. The different appearance of 



the A propagators in (1) and (2) now leads to very differ- 
ent values. In the case of the conventional ansatz, eq. (1), 
we get i?(1440) = 3.9(3) which is well within the range of 
the PDG values of 4(2) Q. In case of eq. (2) the differ- 
ence in the energy dependence of the two decay routes is 
much smaller and we obtain only i?(1440) = 1.3(2). We 
note in passing that due to the strong energy dependence 
of this ratio also the appropriate pole position is crucial. 
If instead of the nominal Breit-Wigner mass pole posi- 
tion we use the pole position evaluated from the speed 
plot of ttN phase shifts, namely M N , = 1371 MeV/c 2 
[HQ], then the values for the ratio decrease to 1-2(1) and 
0.7(1), respectively, using eqs. (1) and (2). 

In summary, the new set of differential data for the 
pp — » ppn + ir~ reaction at T p = 775 MeV supports the 
conclusion that this reaction is dominated by the excita- 
tion of the Roper resonance and its decay into the Ntttt 
channels as derived recently [5j from the analysis of the 
first exclusive measurement at T p = 750 MeV. The new 
data set gives a first experimental evidence for a differ- 
ent energy dependence of the decay routes N* — > A7r and 
N* — » Na. The decay branching of N* — » A7r increases 
by 50% relative to that of N* — » N<r, when increasing 
the incident proton energy from T p = 750 MeV to 775 
MeV, or equivalently when increasing the effective Ntttt 
mass from Mm* = 1264 to M n* = 1272. In this very low- 
energy tail of the Roper resonance we find the N* — > Na 
decay to be clearly dominant with i?(1264) = 0.04 and 
i?(1272) = 0.06. These results are independent of the 
ansatz used for the reaction amplitude. Though we ob- 
serve the low-energy region to be represented very well by 
a Na partition — as suggested e.g. in ref. @] — we also 
see a small but rapidly increasing influence of the A-7T par- 
tition. Due to its k\-k,2 dependence the N* — ► A7T route is 
even likely to finally take over at higher energies. Extrap- 
olating to the resonance pole we obtain i?(1440) = 3.4(3) 
and 1.1(2) using ansatz (1) and (2), respectively. Clearly 
this extrapolation is strongly model-dependent. How- 
ever, as we have demonstrated, the pp — > N Ntttt reaction 
offers the opportunity to experimentally map out the en- 
ergy dependence of the N* — > Ntttt decay systematically 
up to the resonance pole by successively increasing the 
incident proton energy — a program which is currently 
pursued at CELSIUS- WASA. 

We acknowledge the continuous help of the TSL/ISV 
personnel and the support by DFG (European Graduate 
School 683) and BMBF (06 TU 987). 



[1] Particle Data Group, Phys. Rev. D66, 1 (2002) 
[2] O. Krehl et al., Phys. Rev. C62, 025207 (2000) 
[3] E. Hernandez, E. Oset and M.J. Vicente Vacas, Phys. 

Rev. C66, 065201 (2002) 
[4] H.P Morsch, P. Zupranski, Phys. Rev. C61, 024002 

(1999) 

[5] W. Brodowski et al., Phys. Rev. Lett. 88, 192301 (2002) 



[6] L. Alvarez-Ruso, E. Oset, E. Hernandez, Nucl. Phys. 

A633, 519 (1998) and priv. comm. 
[7] H. Calen et al, Nucl. Instr. Meth. A379, 57 (1996) 
[8] J. Patzold, doctoral thesis, Universitat Tubingen 2002; 

http://w210.ub.uni-tuebingen.de/dbt/volltexte/2002/550/ 
[9] R. Arndt et al., Phys. Rev. C56, 635 (1997); 

program package SAID 



6 



[10] F.H. Cverna et al., Phys. Rev. C23, 1698 (1981) [14] D.C. Brunt et al., Phys. Rev. 187, 1856 (1969) 

[11] D.R.F. Cochran et al., Phys. Rev. D6, 3085 (1972) [15] J. Johanson et al., Nucl. Phys. A712, 75 (2002) 
[12] F. Shimizu et al., Nucl. Phys. A386, 571 (1982) 
[13] L.G. Dakhno et al., Sov. J. Nucl. Phys. 37, 540 (1983) 



